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Abstract 

Leptoquark-Higgs interactions induce mixing between leptoquark 
states with different chiralities once the electro- weak symmetry is bro- 
ken. In such LQ models Majorana neutrino masses are generated at 
1-loop order. Here we calculate the neutrino mass matrix and explore 
the constraints on the parameter space enforced by the assumption 
that LQ-loops explain current neutrino oscillation data. LQs will be 
produced at the LHC, if their masses are at or below the TeV scale. 
Since the fermionic decays of LQs are governed by the same Yukawa 
couplings, which are responsible for the non-trivial neutrino mass ma- 
trix, several decay branching ratios of LQ states can be predicted 
from measured neutrino data. Especially interesting is that large lep- 
ton flavour violating rates in muon and tau final states are expected. 
In addition, the model predicts that, if kinematically possible, heavier 
LQs decay into lighter ones plus either a standard model Higgs boson 
or a Z°/W ziz gauge boson. Thus, experiments at the LHC might be 
able to exclude the LQ mechanism as explanation of neutrino data. 



1 Introduction 



Leptoquarks (LQs) appear in many extensions of the standard model. First 
discussed in the classic papers by Pati & Salam |lj and Georgi & Glashow [2] , 
LQs are a common ingredient to grand unified theories [3]. □ They can also 
appear in composite |5J as well as in technicolour models P, [7] . Also in super- 
symmetric models with R-parity violation scalar quarks have leptoquark-like 
interactions [8]. From a low-energy point of view, however, LQs are best 
described in a "model-independent" way, using a LQ Lagrangian based only 
on the minimal assumptions of (a) renormalizability and (b) standard model 
(SM) gauge invariance [9]. An exhaustive list of limits on such LQs from 
low-energy experiments can be found, for example, in |10j . 

Direct searches so far have not turned up any evidence for LQs The 
best limits on pair produced LQs currently come from the DO [12] and CDF 
[T3] experiments at the Tevatron. These typically give limits on LQ masses 
in the ballpark of tulq > (200 — 250) GeV, depending mainly on the final 
state decay branching ratios and on the lepton-quark generation, to which 
the LQ state couples. Considerably more stringent limits are expected from 
the LHC experiments. Depending on the accumulated luminosity, the LHC 
should be able to find LQs up to masses of order of ttilq ~ (1.2 — 1.5) TeV 

Solar [15], atmospheric [TB] and reactor [T7| neutrino oscillation experi- 
ments have firmly established that neutrinos have mass and non-trivial mix- 
ing between different generations. In the SM neutrinos are massless. How- 
ever, non-zero neutrino masses can easily be generated and the literature is 
abound in neutrino mass models [18J. Certainly the most popular way to 
generate neutrino masses is the seesaw mechanism [T$l [2111 |2"TI |2"2"] , countless 
variants exist. However, it is also conceivable that the scale of lepton number 
violation is near - or at - the electro-weak scale. To mention a few examples, 
there are supersymmetric models with violation of R-parity [81 [23] , models 
with Higgs triplets [22J or a combination of both [24J . Also purely radiative 
models have been discussed in the literature, both with neutrino masses at 
1-loop [2S1 EE] or at 2-loop [23 [2E1 [2S1 ED] order. Radiative mechanisms 
might be considered especially appealing, since they generate small neutrino 
masses automatically, essentially due to loop suppression factors. 

In this paper, we study the generation of neutrino masses due to loops 

1 For a recent example of a GUT model with light leptoquarks, see [4]. 
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involving light leptoquarks, in a model with non-zero leptoquark-Higgs inter- 
actions [31J. LQ-Higgs interactions lead to mixing between LQs of different 
chiralities (and lepton number) once electro-weak symmetry is broken and 
thus can contribute non-trivially to the Majorana neutrino mass matrix at 
1-loop level. As discussed below, the peculiar structure of leptonic mixing, 
observed in neutrino oscillation experiments, enforces a number of constraints 
on the LQ parameter space. The main result of our current work is, that 
these constraints can be used to make definite predictions for different decay 
branching ratios of several LQ states. Therefore, the hypothesis that LQ 
loops are responsible for the generation of neutrino mass is testable at the 
LHC, if LQs have masses of the order of 0(1) TeV. 

Before proceeding, a few more comments on LQs might be in order. First, 
for the LQ model to be able to explain neutrino data, non-zero LQ-Higgs 
interactions are essential. Limits on these couplings, on the other hand, 
can be derived from low-energy data such as, for example, pion decay [31]. 
Especially stringent are limits from neutrinoless double beta decay [32J and 
from the decay K° — > e^fi^ [33] ■ However, as we will discuss below, the small 
neutrino masses themselves are up to now the most sensitive low-energy probe 
of LQ-Higgs mixing terms. 

Second, it should be mentioned that LQ loops as a source of neutrino 
mass have been discussed previously in [31] . We will improve upon this work 
in several aspects: (i) We will present neutrino mass formulas containing all 
possible LQ loops, while in [31] only down-type quark loops were considered, 
(ii) [31] concentrated on upper limits on LQ parameters, which can poten- 
tially be derived from observed neutrino masses. We, on the other hand, 
identify the regions of LQ parameters were the neutrino mass matrix is dom- 
inated by LQ loops, thus providing a potential explanation of oscillation 
data. And, lastly but most importantly, (iii) we discuss possible accelerator 
tests of the LQ hypothesis of neutrino masses, to the best of our knowledge 
for the first time in the literature. 

Finally, it should also be mentioned that LQs can be either scalar or 
vector particles. We consider only scalars in details. However, we note that 
most of our results straightforwardly apply also for vector LQs. 

The rest of this paper is organized as follows. In section (T5]) we define the 

2 This is not an entirely new subject. Majorana neutrino masses due to loops involving 
coloured scalars (and vectors) have been discussed first in [27]. However, because in [27] 
it was assumed that these particles have masses M > 10 8 GeV, their contribution to M. v 
was deemed negligible. 
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leptoquark interactions, both with quark-lepton pairs and with the SM Higgs 
boson and discuss the LQ mass matrices. In section fl3]) we calculate the 1- 
loop neutrino mass matrix in the LQ model. Some particularly simple and 
interesting limits are defined and discussed analytically. The typical ranges 
of LQ parameters, required to explain current neutrino data, are explored. 
We then turn to the phenomenology of LQs at future accelerators in section 
(j3J). It is found that some fermionic LQ decays trace the measured neutrino 
angles and thus can serve, in principle, as a test of the LQ model. Next we 
discuss LQ decays to the SM Higgs and to gauge bosons. Higgs (and Z°) 
decays should occur, if kinematically possible, due to the non-zero LQ mixing 
required to explain neutrino masses and thus form a particularly interesting 
signal of the LQ model. We then close the paper with a short summary. 



2 Leptoquark basics 

2.1 Scalar leptoquark Lagrangian 

The SM symmetries allow five scalar LQs. Table Q] shows their SU(3) C x 
SU(2)l x U(1)y quantum numbers, as well as their standard baryon and 
lepton number assignments. LQs which couple non-chirally are strongly con- 
strained by low-energy data [10]. Thus, the states Sq and Sq (as well as S^ 2 
and S^, 2 ), which have the same SM quantum numbers, but couple to (quark) 
doublets and singlets, respectively, are usually assumed to be independent 
particles. Under these assumptions, the most general Yukawa interactions 
(LQ-lepton-quark) induced by the new scalar fields are given by |9] 

+Xf /2 dP L l S{ /2 + Ag } fP L iT 2 l Sf f + Ag) 2 qPmr 2 e Sf/ 2 
+Ag ) ^P L zr 2 5j/ + h.c. (1) 

Here we used the conventions of [10J. Note that eq. is written in one- 
generation notation. In general, all A's are 3x3 matrices in generation space. 
q and I (u, d and e) are the quark and lepton SM doublets (singlets), are 
the scalar LQs with the weak isospin % = 0,1/2,1 coupled to left-handed 
(j = L) or right-handed (j = R) quarks respectively. Thus, in total eq. (pQ) 
contains seven LQ fields. We have also defined Si = r ■ Si. 
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Table 1: Standard model quantum numbers of the scalar leptoquarks. The 
indices 0, 1/2, 1 indicate the weak isospin. The weak hypercharge is normal- 
ized according to Y = 2(Q em — T 3 ). 

The most general renormalizable and gauge invariant scalar LQ interac- 
tions with the SM Higgs doublet (H) are described by the scalar potential 

m 

V = hf Hir 2 A/2 Si + h Sl HiT 2 & Si/2 + Y^ /2 (Hir 2 S\ j2 ) (s\ /2 H) 
+Y Sl (Hit 2 S\h) So + 4° (^i^) 

-(Ml - gl li2) H ] H)^^ + h.c. (2) 

Here is a cumulative notation for all scalar LQ fields with i = L,R (the 
same for z lj2 ) • The diagonal mass terms M|$^$ can be generated by spon- 
taneous breaking of the fundamental underlying symmetry down to the elec- 
troweak gauge group at some high-energy scale. The subsequent electroweak 
symmetry breaking produces non-diagonal LQ mass terms which, in addition 
to the diagonal terms given in eq. (j5J), define the LQ squared- mass matrices. 
These will be discussed next. 

It is important to note, that the first two terms of the scalar potential 
in eq. (J2J) violate total lepton number by two units AL = 2 and, therefore, 
generate Majorana neutrino masses after electro-weak symmetry breaking. □ 
In the limit where and hs t vanish, neutrino masses vanish as well. 

3 Note that an extension of the scalar potential in eq. ([2]) to include the LQ trilinear 
self-interaction terms has been discused in [35] . Such terms, however, introduce violation 
of baryon number and, hence, proton decay, for details see [35] . Therefore we will not 
consider LQ self-interaction terms here. 
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2.2 Scalar leptoquark mass spectrum 



There are four squared-mass matrices which determine the masses of LQs 
with the same electric charge (Q = —1/3,-2/3,-4/3,-5/3). In the in- 
teraction eigenstate basis, defined by S-i/a = (Sq , Sq , S^ 2 , Si), S , _ 2 /3 = 

(Si/2,Sy 2 ,Sy 2 ,S{), 5-4/3 = (So, Si) and 5_ 5 / 3 = (Sy 2 ,S^ /2 ), the squared- 
mass matrices read 



Mia 



h L So v 
Ms 

•31/2 



4 L V\ 

h Sl v 



(3) 



and 



(M 



•2/3 



•S'l/2 



m: 



-4/3 



■5/3 



Y<£ t; 2 

•31/2 



J l/2 



M 



So 



Mat 



1/2 



•31/2 

(LR) 2 

M cR 



V2h Sl v\ 



M\ J 




(LR) V 2' 



(4) 



(5) 



(6) 



1/2 



Here M 



M| — g$v 2 and only the elements above the diagonal have 



been written since the matrices are symmetric, v is the SM Higgs vacuum 
expectation value, v 2 = (2y2G f)~ j • The mass eigenstate basis is defined as 



(2V2G F) 
(S Q )i 



(7) 



where BP is a rotation matrix. The diagonal squared-mass matrices are 
found in the usual way: 



(M%)^ = E?MU&f. 



(8) 



Phenomenological implications of the LQ interactions given in eqs. ([I]) and 
(J2J) have to be derived in terms of the mass eigenstates. For the LQs with 
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charge Q = —4/3 and Q = —5/3, simple analytical expressions for the 
eigenvalues and rotation angle can be found. These are given by 



Ml, = \ (m 2 xx + M\ 2 - y^ + ^-Mf^) , (9) 



and 



2M 2 12 



tan2fl 12 = 2 " . (10) 



Here, M xl , M| 2 and Mf 2 stand symbolically for the corresponding entries in 
the mass matrices eqs (jSJ) and For LQs of charge Q = —1/3, —2/3 we 
will diagonalize the mass matrices numerically below. However, the following 
approximate expressions are useful for an analytical estimate of parameters. 
The rotation matrices which relate the interaction and mass eigenstates can 
be parametrized by six rotation angles, namely, 

RQ = R Q {e M )R Q (9 M )R Q (6 14 )R Q {6 23 )R Q {6 13 )R Q (e 12 ). (11) 

In the limit where the off-diagonal entries in the mass matrices eqs ([3]) and 
(j4j) are smaller than the difference between the corresponding diagonal ones, 
it is possible to find approximate expressions for the rotation angles also 
in this more complicated case. As discussed in the next section, for the 
neutrino masses the most relevant angles are 6^ -2 ^ 3 , 6^ -1 ^ 3 and #^ _1 ^ 3 . For 
the angles in the Q — 1/3 case on can use eq. ffTUj) as an estimate, with 
obvious replacements of indices. For the angle 6^ -2 ^ 3 , however, since the 
relevant M| 4 = in the mass basis, a more complicated expression results: 



=2/3 ^Y S R lf M^ 



(M s n 2 -M §i/2 )(M Sl -^ s n 2j 



Eq. dT2]) is exact in the limit Y£ 2 = gf R) =0. It remains a reasonable 



c {LR) 

1/2 ' " 

LQ mixing angles are small numbers. 



(factor-of-two) estimate as long as .g^^ — ^s^^i A*' an d a ^ Q = 



6 



(£-2/3)5 (^2/3) 



(5-1/3)] 



(a) 



(6) 



i/i/ 



Figure 1: Feynman diagrams for Majorana neutrino masses. Diagram (a) 
[(b)] give contributions to the neutrino mass matrix from it-type [d-type] 
quarks loops. 

3 Neutrino masses from Leptoquark loops 
3.1 Analytical formulas 

LQ-qv Yukawa interactions can be derived directly from the Lagrangian ([1]). 
In the interaction eigenstate basis, they have the following form 



Clq^ = Xl /2 uP L u(S^_ 2/3 + \%u*P L v{S$_ m + h.c. (13) 



and 



C'LQ-du = Xf dP L v(S^)l 1/s ~ \ L So dcP L v(S5)\ /3 + \%#PlM(Si% + h.C. 

(14) 

Rotating to the mass eigenstate basis, the non-trivial mixing among LQs 
from different SU(2)l multiplets lead to neutrino Majorana masses at one- 
loop order as shown in figure [TJ A straightforward calculation of the Majo- 
rana neutrino mass matrix from these diagrams gives 



M u = M7 + M 



down 



(15) 



where the matrix Ai^ p from diagram (a) reads 

3 



(M 



v Ju' 



— 2 £ m k B (0, mlM])Rf R f 



J=l-4 
k=u,c,t 



X 



(^Su^iki^sji'k + {^S 1/2 )i'k{X Sl )ik 



(16) 



Here R 2 ^ 3 is the rotation matrix that diagonalizes the mass matrix of Q = 
—2/3 LQs, eq. fll]), and B o (0, m 2 ,, M 2 ) is a Passarino-Veltman function 
The matrix J^i^ own from diagram (b) can be written as 

[K-) lV = ^ £ m k B (0,mlM 2 )R% 3 



3=1-4 
k=d,s,b 



x 4 i?^ 3 



(A i 5 i/2 )ifc(A 5i ) i / fc + (A« 2 )i'fe(A Sl )ife 
+ i? J 1 { 3 (A| ' )ik(^s )i'k + (A| 1 )i'fe(As )ifc }• (17) 



liJ 1/2 )ifc( A i )i'fc + (Af i/2 ) - ' xL 

Here R 1 ^ 3 is the rotation matrix that diagonalizes the mass matrix of the 
Q = — 1/3 LQs given in eq. ([3]). Note, that in the limit of unmixed LQs, i.e. 
Rij = 5ij, the neutrino mass matrix vanishes. 

The Passarino-Veltman function Bq [36] contains a finite and an infinite 
part. However, since the LQ model does not have a neutrino mass at tree- level 
there are no counter terms, which allow to absorb infinities. The infinite parts 
of the B functions therefore must cancel among the different contributions 
in eq. (fTBl and eq. (fT7l) . Using the parameterization of the LQ rotation 
matrices given in eq. ( TTT1) . we have checked that this is indeed the case. The 
resulting formula can be expressed as a sum of differences of Bq functions 
only, thus cancelling all infinities. Since the coefficients in these formulas are 
rather lengthy (and of little use), we will not give them explicitly. 
Diagonalizing eq. ffl5l) gives the neutrino masses and mixing angles, 

U T M U U = Mf ag . (18) 

In standard parameterization U is written as 

s 13 e~ iS 

U = I c 23 s 23 I I 1 M -sia c 12 I , (19) 





where c^- = cos% and = sin 6*^ and 5 is a CP- violating phase. Since 
we will consider only real parameters, 5 = 0, n and we have not written any 
Majorana phases in eq. fTT9|) . 

In general, the neutrino mass matrix receives contributions from diagrams 
involving up-type (it-loops) and down-type (d-loops) quarks. In order to find 
the eigensystem of eq. ( Tl5i) one has to solve a cubic equation. However, much 
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simpler analytical formulas can be derived, if one particular loop dominates 
over all others. For example, in the limit where only the top loop contributes 
to A4 U one finds Det [M v ] = 0, i.e. one of the three eigenvalues of the mass 
matrix goes to zero. Note, that in this limit the model can produce only a 
(normal) hierarchical neutrino spectrum. 

Analytical expressions for the two non-zero neutrino masses can be found 
easily in the limit Det [M v ] =0. It is useful to define two vectors in param- 
eter space 

R= [(A?)i,(A?MA?) 8 ], 

L= [(A£)i,(A£MA£) 3 ]. (20) 

Here, (Af' j; )j = (Xfp)jk, with j being the leptonic index, whereas we have 
suppressed for brevity the hadronic index k. The indices / and I' stand 
symbolically for I = Si/ 2 and I' = Si if the top loop dominates, or I = S%/2 
and I' = Si or I' = S , if one of the bottom loops dominates. In terms of 
these vectors the two non-zero neutrino masses are given by 

=.F(|R.L|=f|R||L|), (21) 

where T is given by 

^-^E m k B o {0, ml M])R%R9, (22) 

71 3=1-4 

Q = 1/3,2/3 and s = 1,4, depending on which contribution to M. v is most 
important. The ratio between the solar and the atmospheric scale is thus 
simply given by 



H (\ 


R L 




R| 


|L| 


Am 2 23 V 


R L 


+ 


R| 


|L| 



^7^^ D t I "d t ) • (23) 



Note that R is independent of J 7 , i.e. independent of LQ masses and mixings. 
Relations among neutrino mixing angles and the Yukawa couplings can be 
found by using the eigenvalue equation for the massless neutrino [291 SO] , 

M u v = (24) 

where the eigenvector t> is given by 

v*= / W) ■ (25) 
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Solving eq. (1211) yields the result 



mi 2 m 33 - m l3 m 23 , m 12 m 23 - m l3 m 22 

e= 2 — , e = — , (26) 

m 22 m 33 - m z 23 rn 22 m 33 - m\ 3 

where are the entries of the neutrino mass matrix M. v . Interestingly, 
eq. (I2"6"j) can be expressed in terms of neutrino angles only. For a normal 
hierarchical spectrum, i.e. m ui23 ~ (0,m, M), where M (m) stands for the 
atmospheric (solar) mass scale one obtains 

e = tan 8 12 C ° S ^ 23 + tan 8 13 sin 8 23 , (27) 

COS #13 

e = tan6>i 2 — — — tan# 13 cos# 23 . (28) 
cos 6> 13 

On the other hand, the expressions for e and e' eq. (126]) depend on the 
entries in the neutrino mass matrix which are determined by the LQ Yukawa 
couplings, 

(Af)3(A£)i-(Af)i(Afr)3 



(A?) 3 (A£) 2 -(A*) 2 (A£) 3 
(Af) 2 (Af>)i-(Af) 1 (Af',) 2 



(29) 



C ~ (Af)3(A^) 2 -(Af) 2 (A^,) 3 - W 

The above equations allow to relate the Yukawa couplings directly to the 
measured neutrino angles. Note also, that current neutrino data require 
both, e and e', to be non-zero. 



3.2 Neutrino data and parameter estimates 

Before discussing the constraints on LQ parameter space imposed by neutrino 
physics, let us briefly recall that from neutrino oscillation experiments two 
neutrino mass squared differences and two neutrino angles are by now known 
quite precisely |37j. These are the atmospheric neutrino mass, Am\ tm = 
(2.0 - 3.2) [10- 3 eV 2 ], and angle, sin 2 9 Atm = (0.34 - 0.68), as well as the 
solar neutrino mass Amq = (7.1 — 8.9) [10~ 5 eV 2 ], and angle, sin 2 # Q = 
(0.24 — 0.40), all numbers at 3 a c.l. For the remaining neutrino angle, the 
so-called Chooz [38] or reactor neutrino angle 9r, a global fit to all neutrino 
data [37] currently gives a limit of sin 2 Or < 0.04 @ 3 a c.l. 

Neutrino oscillation experiments have no sensitivity on the absolute scale 
of neutrino masses, but the atmospheric data requires that at least one neu- 
trino has a mass larger than M = m^ tm > 50 meV. The minimal size of LQ 
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Yukawa couplings and LQ-mixing, required to explain such a neutrino mass, 
can be estimated from eq. (1211) . Parameterizing the rotation matrices as in 
eq. ffTTj) we can estimate T as 

T ~ -^— 2 m k sm(20 3s )AB 3s . (31) 

l07T Z 

Here, sin(2#q s ) stands symbolically for the largest LQ mixing angle, and 
AbJ = 1(0, Ll m!) - Zo, ml k). The finife J*B> of B \ is fiven b y 

Q // n 2 2n ^ l og (m2) - m 2 log(m 2 ) 

" L k " b j 

The maximum allowed value of \ABij\ for tulq < 1-5 TeV is |A5y| ~ 3 (3.5) 
for rrik = m t (m& = 0). With the current central values for the quark masses 
[39] we then find the maximal value (s) of T as 

jrmax _ [ 49 . Q14 . 41 . 10 -2 . 3. 10 -3 . 1.7 . 10~ 4 : 8-10~ 5 ] GeV, 

for t : b : c : s : d : u (33) 

for maximal LQ mixing, i.e. sin(26 , 3s ) = 1/2. For this value of T the min- 
imum values for the Yukawa couplings required to explain the atmospheric 
mass scale are then very roughly given by 

(AflUMA^Vt > 5-1 x 1(T 12 , (Af i/2 WA£ 01 U > 1-8 x 10- 10 (34) 
(Af 1/2 )ic(A£> c > 6.0 x 10- 10 , (Xl /2 )is(X L s Ji's > 8.0 x 10" 9 
(Af i/2 ) id (A| i )^ > 1.5 x 10- 7 , (Af i/2 ) m (A§>„ > 3.0 x lO" 7 

Obviously, unless the \n~ follow a hierarchy inversely proportional to the SM 
quark masses, third generation quark loops give the by far largest contribu- 
tion to the neutrino mass matrix. 

We should compare the minimal values of eq. (134")) with the constraints 
coming from low-energy phenomenology. The most stringent upper bounds 
for the first generation Yukawa couplings ((Af i/2 )ii and (A^Ja) are currently 
found from the upper limit on the lepton flavour violating process \i Ti — > e Ti 
MM 

^ ^ ^ < 2.6 xHT 7 f M l V 



( A 5 1/2 )ll( A 5 1/2 )21 

(A^)n(A^) 21 < 1.7 x 10- 7 ( 77 ^ ZT7 ) • (35) 



ViooGevy ' 

(-^-) 

V lOOGeV / 
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Here j labels the corresponding mass of the LQ eigenstate. Upper bounds 
for the second (and third) quark generation couplings come from the charged 
lepton flavour violating decay \i —>■ cy and are given by 



(A3 1/2 Jl2(13) 1/2 j22(23) < 



(A S J 12(13) (A Sl ) 22(23) < 



1.8 X 10 



1.8 x 10 



,-5 



-5 




Here, we have updated [10] with the current experimental upper limit on 
Br(// -+erj) [39]. 

Although eq. ()35p constrains a different combination of left- and right- 
LQ couplings than eq. (1341) we conclude that, barring cases where some fine- 
tuned cancellation between different LQ contributions occur, we expect that 
first generation quark loops can not explain current neutrino data. Second 
and third generation LQ loops, on the other hand, could both produce the 
observed neutrino masses, consistent with all phenomenological constraints. 
However, considering the hierarchy in m c /m t ~ 8 • 10~ 3 and m s jm^ ~ 0.02, 
from now on we will concentrate on third (quark) generation LQs. Note that, 
comparing eq. (l34j) with eq. (!36|) one finds that the atmospheric mass scale 
can be generated consistent with low-energy constraints for LQ mixing as 
small as 10 -6 (10 -5 ) in case of top-loops (bottom- loops). These numbers are 
significantly smaller than constraints derived from other low-energy processes 



The observed large mixing angles in the neutrino sector require certain 
ratios of Yukawa couplings to be non-zero. This can be most easily under- 
stood as follows. One can use eqs ([TBI and (Tl9l to invert the problem and 
calculate the neutrino mass matrix in the "flavour basis" (in the basis where 
the charged lepton mass matrix is diagonal). The resulting }A V in the general 
case is a complicated function of the eigenvalues and mixing angles. However, 
as first observed in [5T], the so-called tri-bimaximal mixing pattern 



is a good first-order approximation to the observed neutrino angles. In case 



PI S3]. 




o \ 



i 



(37) 
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of hierarchical neutrinos j\A d J a9 = (Q,m,M) it leads to 

M* PS = - I M -M \ + - I m m m I . (38) 




Comparing eq. (1581 with the index structure of eq. (TIBl) and (flTl) . one expects 
that 

(AfMAfk + MMA?)! « (A?) 2 (Af) 8 + (A? ) s (Af ) 2 (39) 
(Af , ) 2 (A? ) 8 + (Af)a(A? ) 2 > (Af) 2 (Af ) 2 - (Af ) 3 (Af ) 3 (40) 

for the couplings which give the largest contribution to M. v . Eq. (l39l is 
essentially due to smallness of the reactor angle, while eq. fHUj) follows from 
the observed near-maximality of the atmospheric angle. Note that, if more 
than one loop contributes to M. v of eq. ( Tl5l) . m vx ^ 0, but the "large" off- 
diagonal entry in the (2,3) element of M. v always requires (Aj^) 2 ~ (X^ R ) 3 , 
for at least one LQ state. Finally, it should also be mentioned that the 
smallness of solar versus atmospheric splitting requires that the vectors R 
and L, defined in eq. fl20l . are nearly aligned for all vectors contributing to 
A4 U , compare eq. ( |23l) . 

There are three different contributions to the neutrino mass matrix, see 
eqs ( 1T6|) and ( ITT]) . The top loop is proportional to O®^ 2 ^ ~ ^^ 1/2 ' hs 1: while 

the bottom loop is either proportional to 6®^ 1 ^ ~ h Sl or to 6*^ -1 ^ 3 ~ hg. 
Lacking a theoretical ansatz for these parameters, it is not possible to predict 
which of these give the dominant contribution to the neutrino mass matrix. 
However, since rrib/mt ~ 2%, the top-loop will be most important, if all LQ 
mixing angles (and Yukawa couplings) are of similar size. We will refer to 
this case, M. v = M.\,, as scenario I. O^ 2 ^ 3 , on the other hand, can be much 
smaller than the corresponding angles in the down-type loops in those parts 
of parameter space where all relevant off-diagonal entries in the LQ mass 
matrices are small. In this case M. v ~ A4 b u and we will refer to this situation 
as scenario II (II. a: if hg h Sl /v and II. b: if hs x /v <C hg ). 



4 Leptoquark collider phenomenology 

LQs, once produced, will decay almost instantanously. There are two differ- 
ent sets of possible final states. In the current model, apart from the usual 
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decays into a quark and a lepton there are also vector (W and Z°) and 
scalar (h°) final states, if kinematically allowed. We will discuss first the 
fermionic decays. 

4.1 Fermionic LQ decays 

Fermionic decays of the LQ mass eigenstates are dictated by the Yukawa 
interactions given in the Lagrangian ([1]). Possible final states can be either 
iiiik, lidk, ViU k or Vid k . Most interesting, from the phenomeno logical point of 
view, are final states with charged leptons, since these allow to tag the flavour. 
Partial widths for two-body final states can be calculated in a straightforward 
manner. For charged lepton final states these are given by, 

r[(3_ 1/3 ), - m = < 1/3 {[(X L Sl hR% 3 } 2 + [(X L s hR]i 3 ] 2 + [(Af ) i3 4 2 /3 ] 2 } ; 

(41) 

r[(S_ 2/3 ), - Ob] = ni 2/3 {[(Af i/2 ) i3 4 /3 ] 2 + m i/2 hR% 3 } 2 } ; (42) 

n@-w)i - m = <4 /3 {t(Af>< 3 ] 2 + i(xihR% 3 ?} ; (43) 

T[(5_ 5/3 ), - m = ni 5/3 { {(X L Sl/2 hR% 3 } 2 + [(Af i/2 )„ 3 4 /3 ] 2 } • ( 44 ) 
Here, Uq is an overall constant given by 

*3 r 771 ^ -\- 771 ^ 1 

^ = t^— 1 - A 1/2 K„ ml m?) (45) 

v IbimiSj L Trig \ 3 

with A 1 / 2 (a, b, c) the usual phase space factor, A(a, b, c) = (a+6— c) 2 — 4a6, and 
m Sj , m q and the corresponding LQ, quark and lepton masses. Absolute 
values for the LQ widths can not be predicted. However, minimal [maximal] 
values can be estimated from the atmospheric neutrino mass scale [low-energy 
bounds]. Putting all parameters to their extreme values fermionic widths 
could be as small [large] as 0(eV) [O(MeV)]. 

In the above equation, we have written only the partial widths to top 
and bottom quarks. Formulas for the lighter generations can be found with 
straightforward replacements of indices. However, since the widths, eqs (1411) 
- ( 1441 . are not suppressed by quark masses, our assumption that 3rd genera- 
tion quark loops give the dominant contribution to Ai v , can in principle be 
checked experimentally. For example if Br(S^ 5 / 3 — > t + £\ k)/B r {S-5/3 ~^ 
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3 + J2i h) and Br{S-±/z ~^b + J2i h)/Br(S-4/3 -> j + £\ k), where j stands 
for any non-b jet, are larger than m c /m t , charm (and up) loops are guaran- 
teed to be sub-dominant. Similar tests can be devised for the case of bottom 
loops. 

Note that if the mixing between different LQs is small, as is generally 
expected, the decays of some of the LQ states Sq are controlled by the same 
Yukawa couplings that determine the non-trivial structure of the neutrino 
mass matrix. This observation forms the basis of the different decay pattern 
predictions discussed below. However, one complication arises from the fact 
that we can not predict if the decays of the lightest or one of the heavier of the 
LQ states is dictated by the Yukawas fixed by neutrino physics. Again, in the 
limit of small LQ mixing, this question can be decided experimentally, in prin- 
ciple. Consider, for example, scenario I, M. v ~ .M* . The decays controlled 
by neutrino physics are those governed by A5 and A^ . These states couple 
mainly to lepton doublets. Their components have the same diagonal entries 
in the LQ mass matrices, we expect them to have similar masses. These 
states should have very roughly r[(5 , _ 5 / 3 ) — > ^2£it\ ~ r[(S , _ 2 /3) — > ^2^it\ 
and r[(£_4/ 3 ) -> ~ r[(£_i /3 ) -> The other Q = 4/3,5/3 

mass eigenstates mainly couple to singlet leptons, i.e. these state do not 
decay to neutrinos. In what follows below, we will always assume that the 
small mixing limit is realized and the LQ states relevant for the experimental 
cross-checks can be identified. 

Combining eqs (1231 and (l4"0l) with the decay rates eqs (jHj) - (|4"4"|) , one can 
derive some qualitative expectations for some ratios of branching ratios of 
fermionic LQ decays. In general, the constraint from the large atmospheric 
angle, plus the smallness of R = Amg/Am^ tm can only be full-filled if there 
are (at least) two LQ states which have similar branching ratios to muonic 
and tau final states. At the same time, these LQ states should have less final 
states with electrons, essentially due to eq. (1391) and the upper limit on the 
reactor angle. 

Much more detailed predictions for fermionic decays of LQs can be made 
in the explicit scenarios defined in the last section. We will first discuss in 
some detail the results for scenario I, M. v = }A l v . For all figures presented 
in the following we have scanned the Yukawa parameter space randomly, 
in such a way that all low-energy bounds are obeyed. We then numerically 
diagonalized the resulting neutrino mass matrices and checked for consistency 
with current neutrino oscillation data. Different correlations among ratios of 
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Figure 2: Ratio of decay branching ratios yj Br^ 3 Br^J 3 / yj B 

tan 2 ^23- Vertical lines indicate current 3a range for tan 2 823 while horizontal 
lines determine the predicted range for this observable. 



branching ratios with the different pieces of neutrino data are then found. 

Figure [2] demonstrates that yj Br^ 3 Br^J 3 / yj Br^j 3 Br^ 3 is correlated with 

the atmospheric mixing angle. For the best fit point value tan 2 623 = 1 one 

expects yj Br^ 3 Br^J 3 ~ yj Br£j 3 Br%j 3 . Using the current 3a range for the 

atmospheric mixing angle this observable can be predicted to lie within the 
interval [0.4,4.7]. 

We have found that there exists an upper bound on the ratio of branching 
ratios ( 

Br f/3 Br l/3 

< 9 x 10~ 2 (46) 



H% Br l% + yjBr% 3 Br% 3 

which can be derived from the ratio R = Am^ / Am\ tm as shown in figure [3j 
This bound shows that the product of branching ratios Br^Br^ is expected 

to be nearly two orders of magnitude smaller than the sum of Br^^Br £, 3 and 

f> r tr r> br 
-° r 5/3 r 4/3- 

Individual values for electron final state decay branching ratios are shown 
in figure HI It can be seen that the smallness of Br^, 3 Br\ e , 3 can be due to 
the smallness of either Br b ^ 3 or Br^ 3 . This implies that for one of the two 
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LQ eigenstates (Q = 4/3 and Q = 5/3) electron final states could be as 
large as ~ 20%, but only if the other LQ state shows a very much supressed 
branching ratio to electrons. 

Numerically we have found that there is certain combination of ratios 
of branching ratios that is correlated with sin Or = sin #13 as shown in fig- 
ure [5J With the current upper limit on sin 9r, this ratio is not very much 
constrained. However, a future measurement of sin#R, would confine this 
ratio to lie in a very small, albeit double-valued, interval and thus such a 
measurement could become a powerful experimental cross-check of the sce- 
nario discussed here. Note also that this ratio approaches 1 for small values 
of sin#R, thus also an improved upper limit on this angle will lead to an 
interesting constraint. 

Finally, from equation (123]) one expects that 

^=, T Brf /3 Brf /3 + sjY~~ T Br% 3 Br% 3 ' 

The neglection of electron final states in the 2nd equation above is motivated 
by eq. (j4"Sj) . Numerical results are shown in figure [6j The spread of the points 
in the plot gives the precision with which the ratio Br-/Br + can be predicted, 
neglecting electron final states and scanning over the allowed ranges of other 
neutrino physics observables. As demonstrated by figure the observable 
Br_/Br + is currently expected to lie within the range [7.5 x 1CT 3 , 2.9 x 10~ 2 ]. 

All results shown in the figures (j2J)- ([6]) are based on the assumption that 
the top-loop gives the most important contribution to the neutrino mass 
matrix. However, very similar results can be obtained if the bottom loop 
dominates in either scenario II. a or scenario II. b. We will not repeat the 
discussion in detail here. The results for scenario II. a can be obtained by the 
replacement of S5/3 — > S2/3 and scenario lib by the replacements S5/3 — >■ S2/3 
and S4/3 — > S1/3 in all equations and figures above. 

In summary, qualitative expectations for some ratios of branching ratios 
of fermionic LQ decays can be derived from the requirement that LQ loops 
explain neutrino oscillation data. In general, lepton flavour violating decays 
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with similar branching ratios to muonic and tau final states are expected for 
some specific LQ decays. Sharp predictions for various decay modes can be 
made, under the reasonable asumption that one LQ loop dominates over all 
others. 



4.2 Leptoquark decays to Higgs and gauge boson final 
states 

Since the current lower limit on the mass of a standard model like Higgs boson 
is m h o > 114.4 GeV (22], one expects that LQs can decay also to standard 
model gauge bosons, W ± and Z°, if the Higgs final state is kinematically 
possible. We will therefore discuss partial decay widths to Higgs, W ± and 
Z° final states jointly in this subsection. 

In the model discussed here, heavier LQs can decay to lighter LQs plus a 
standard model Higgs boson, i.e. (Sq)j h° + (Sq)^, due to the interactions 
given in eq. (J2J). Partial decay widths can be written as 

r[{S Q )i - h° + ($,),] = ^ Q m Sj \ l /\l,r ih r h ). (48) 



Tij — 



Here, the arguments of A 1 / 2 (a, 6, c) have been defined dimensionless, 
m|./m|. and = mf i0 /m 2 s .. The effective couplings ]}q for the different 
values of Q = —1/3, —2/3, —4/3, —5/3 are defined as 

(LR) h {L) (L) 

~ _ 9 So v pl/3 p l/3 | n Sp pl/3 p l/3 | K S V pl/3 R l/3 

9 - 1/3 ~ ~W s Ril Ra + 7^s. jl 13 ~^W Rjl R * 4 

, %o pl/3 p l/3 , K S V pl/3 D l/3 , n Sl pl/3pl/3/, Q N 
+ —7= K j2 K i3 + _ K j2 K iA + —fK K j3 K iA I 49 J 

V2m s . 2 m Sj V2m g , 



9-2/3 



yL yR . 

S1/2 V 2/3 p2/3 . S 1/2 V 2/3 p2/3 , I* Si p2/3 p2/3 



(LR) 

9 Sl/2 V 2/3 D 2/3 

+ 2 m s j2 a ' 



p^/o p^/o i -Ji p^/^ p^ 

„ tt n -Kj 3 H -ft,-i -Kj 4 

2 m 5 . 31 m Sl Jl 



9-4/3 



Ysi v 



s/2 m Sj ' 



p4/3p4/3 



(50) 
(51) 
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(LR) 

~ g S 1/2 V D 5/3 D 5/3 

9-5/3 = —Z Rj{ Ri2 ■ (52) 

2 m Sj J 

are the rotation matrices, which diagonalize the LQ mass matrices. Note, 
that the above couplings contain the same parameters which induce neutrino 
masses due to LQ mixing. 

For any given set of LQs of charge Q the couplings with the Z° can be 
written as 

-^-^E ( T 3 - Q sin2 M Si^Syt. (53) 

COS V\y '—~ 

Non-diagonal couplings of the Z° gauge boson to different LQ states of the 
same Q, but different T 3 appear, after rotation to the mass eigenstate basis. 
The partial decay width can be written as 

r l{ S Q)j Z° + (S Qh] = ^(l, r w , r z ), (54) 



where rz = (m^o/mj) 2 and 



O-i/s = ~RfR]l z (55) 

a _ I c.2/3 c.2/3 , 3 p2/3 p2/3x 

0-2/3 - K il + 2 H j* K i± > 

a _ p4/3 p 4/3 

P_4/3 — — -n-j2 -^2 ■ 

Note that Q = —5/3 LQs do not have any decays to Z° bosons, since their 
couplings to Z° are completely diagonal. Closer inspection of eq. (!5"5i) re- 
veals that the decays to Z° states can occur only if LQ mixing (by the same 
parameters which govern the Higgs final states) is non-zero. Thus, also ob- 
servation of Z° final states gives valuable information about the parameters 
ineq. ©. 

Heavier LQs can decay to a lighter one and a W ± gauge boson, Sq — ► 
W + Sq/, where Sq and Sq> are members of the same doublet (triplet). 
Possible decays therefore are: 



(S- 5/3 )i + (5 , _ a/3 ) i , (56) 

(5_ 2/3 ) . <_> W~ + (S_ 1/8 )J, (57) 

(5_ 4 /3), <-> + (5_ 1/3 )<, (58) 
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where the processes in (1561) . (1571) and (15"8j) come from the decays of the mem- 
bers of the doublet 5*1/2, 5*1/2 and the triplet S±, respectively, after rotation 
to the mass eigenstate basis. Note that the process in eq. (|5"T|) can also come 
from the decay of the T3 = 1 to the T3 = components of the triplet. The 
decay widths for the processes in (1561) . (157|) and (1581) can be written as 

n@ Q )i - W ± + (S Q/ U = ^^^(l.ry.Tv). (59) 
Here rw = / m 2 s , and the mixing factors are given by 

0_ 5/3 = (R 2 ^) l2 (R 5 / 3 ) n + (R 2/3 ) l3 (R 5/3 U (60) 
O-2/3 = (R l/3 UR 2/3 ) n + s/2(R 1/3 UR 2/ %4 (61) 
^4/3 = V2(i? 1/3 ) i4 ( J R 4/3 )i2. (62) 

Our formula eq. (1391) agrees with the one calculated earlier in [12], once LQ 
mixing is properly taken into account. 

We now turn to a discussion of typical ranges for the branching ratios 
of bosonic final states. We will first discuss the example of decays of LQs 
with Q = 4/3, assuming the decay to Q = 1/3 LQs plus W ± is kinematically 
closed. 

Fig. ([7]) shows a set of numerical examples of branching ratios of the 
heavier of the Q = 4/3 LQ mass eigenstate to fermionic, h° and Z° final 
states, for some typical choices of parameters, see figure caption. For Yukawa 
couplings of the order A ~ 1CT 3 values of Ys 1 as small as Y$ 1 — 1CT 2 can lead 
to observable branching ratios into bosonic final states. 

While we can not predict whether fermionic or bosonic final states will 
dominate, it is interesting to note that the current LQ model makes a definite 
prediction for the ratio of branching ratios of h° and Z° final states, if m^ -4,73 
is sufficiently larger than m^ -4 ^ 3 + m^o . This can be understood as follows. If 
m Q-4/3 . g muc j 1 i ar g er than m^ _4//3 , m h o and m z o, one can neglect the phase 
space correction factors, \(l,x,y), and approximate the Q = 4/3 mixing 
angle by 6q=a/3 — v / 2Xs 1 t> 2 /(m 2 5_4//3 ) 2 . The ratio of the partial widths to 
Higgs and Z° states is then simply given by 

n$^z°H§*h)„£*„ 8j (63) 



r[(5 4 / 3 ) 2 ^ ^ + (54/3)1] c^mlo 
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Figure 7: Typical values for decay branching ratios for the heavier Q = 4/3 
LQ, in case W ± final states are kinematically closed. Branching ratios are 
plotted versus the "average" Yukawa coupling A = \/$^ A^ 3 , for different 
values of and m^ _4//3 . Full line (red): fermionic final states, dashed line 
(blue) Higgs final state, dot-dashed (magenta) Z° final state. In all figures 
rnQ-4/3 has been set to mf _4//3 = 250 GeV and we have chosen m^o = 115 
GeV, motivated by the LEP limit. Top left: (Y Sl = 0.01,mJ =4/3 = 400 
GeV); top right: (Y Sl = 0.1,mJ =4/3 = 400 GeV); bottom left: (Y Sl = 
0.01,m? =4/3 = 800 GeV); bottom right: (Y Sl = 0.1,m? =4/3 = 800 GeV). 



independent of all non-SM parameters. This explains the ratio observed in 
the numerical examples of fig. (j7|) and constitutes a nice consistency test for 
the LQ model of neutrino masses. 

We now turn to W ± final states. In general, in electro-weak symmetry 
breaking new mass terms for LQs, which are members of the same multiplet, 
could be generated by some non-SM scalars, potentially introducing large 
splitting within a given multiplet. In this case, LQ decays to W ± states 
could occur independent of LQ mixing between different multiplets. Then, 
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since LQ-H^ decays are of order g 2 they would easily become dominant 
once kinematically allowed. In the current model, however, mass splitting of 
LQs within the same multiplet comes only from LQ mixing, see eqs ([3])- ([6]). 
Thus LQ-VV^ final states should have widths similar to the Z° final states 
discussed above. Consider, for example the decays of a Q = 5/3 LQ. The 
mass matrix of the Q = 5/3 LQs, see eq. (jSJ), contains the same parameters 
as a (2-by-2) submatrix of the Q = 2/3 mass matrix, compare to eq. 0. If 
the other Q = 2/3 states are heavier than these states, we can give a similar 
estimate of the ratio of Higgs and W ± final states, as has been derived above 
for Z° final states, see eq. (1631) . Assuming again m^ _5//3 being much heavier 
than final state particles, we find 

r[(g 5/3 ) 2 ^fe° + (g 5 /3)i] ^ L<^ 06 3 (64 ) 
r[(5 B /s) 2 -w± + (5 2/s )i] W v 2 

In the general situation, however, when all Q = 2/3 states are relatively 
light, the branching ratio to h° final states can not be predicted accurately. 
Thus, W final states can not provide an as valuable test for the model as 
is the case for Z° decays. 

In summary, heavier LQs will decay to bosonic final states, if kinemat- 
ically allowed. Since in the current model all these decays are induced by 
the presence of the LQ-Higgs interaction parameters, observing such decays 
are an essential test of the LQ model of neutrino mass. Branching ratios for 
bosonic final states typically fall into the range C(10 -4 — 1), for LQ-Higgs 
couplings order O(10~ 2 — 1). Although we have discussed only the cases 
Q = 4/3 and Q = 5/3, bosonic widths of LQs with other electric charges are 
expected to show a very similar parameter dependence (and therefore similar 
branching ratios). 



5 Summary 

LQ fields with baryon number conserving Yukawa interactions can have 
masses at or near the electro-weak scale. If these LQ fields couple to the 
SM Higgs, the resulting model generates neutrino masses at the one-loop 
level. In this work we have explored the phenomenological consequences 
of LQs as the origin of the observed neutrino masses for future accelerator 
experiments, such as the LHC. 
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Fermionic decays of (some of the) LQ states trace the neutrino angles, 
i.e. certain ratios of decay branching ratios can be predicted from current 
neutrino data. In general one expects that those LQs, which give the dom- 
inant contribution to the neutrino mass matrix, if (pair) produced at the 
LHC decay with sizeable flavour violation. For these states there should be 
a similar number of events with r ± yU =F final states, as there are final states 
with muon and tau pairs. One also expects a smaller number of events of the 
type e ± fi T (and e ± r =F ), although the details in this case are more involved, 
as discussed above. 

In this context we would like to stress that one of the basic assumptions 
applied in practically all accelerator searches for LQs is that LQs couple 
only to one generation of leptons and quarks at a time. As discussed at 
length above, such completely generation diagonal couplings would exclude 
LQ-loops as an explanation of neutrino oscillation data. Extending the LQ 
search to lepton flavour violating decays thus should be considered seriously 
by experimentalists. 

We have also discussed how, in some specific scenarios, much more de- 
tailed predictions can be made. Given the observed hierarchy of standard 
model quark masses, it seems reasonable to assume that contributuins from 
3rd generation quark loops dominate the neutrino mass matrix. For the case 
of top quark dominance, our results are summarized in the figures (|2l)- ([6]). 
Similar results hold in case of pure bottom-loop dominance. 

Finally, an important test of the hypothesis that LQs can generate Ma- 
jorana neutrino masses, is the search for decays of heavier LQs into lighter 
ones plus a standard model Higgs or gauge boson. Any observation of a 
non-zero branching ratio for the decay Si —>■ Sj + h°/Z° constitutes proof for 
LQ mixing, which is the basic ingredient for the LQ explanation of neutrino 
masses. If LQs are found at the LHC, the search for such decays should be 
made a priority. 
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